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Abstract — Sequences sets with low aperiodic auto- and cross- 
correlations play an important role in many applications like 
communications, radar and other active sensing applications. The 
use of antipodal sequences reduces hardware requirements while 
increases the difficult of the task of signal design. In this paper we 
present a method for the computation of the Integrated Sidelobe 
Level (ISL), and we use it to calculate the asymptotic expression 
for the ISL of a set of sequences formed by different rotations 
of a Legendre sequence. 

Index Terms — Integrated Sidelobe Level, antipodal sequences, 
Legendre Sequences, auto-correlation, cross-correlation. 



To facilitate the discussion, denote the sum of squares 
corresponding to the auto-correlation terms as 
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and the sum of squares corresponding to the cross-correlation 
terms as 



I. Introduction 

THE design of sequences set with good correlation proper- 
ties is present in many fields of engineering such as radar, 
sonar, communications, medical imaging and so on. Good 
auto-correlation properties means that any sequence in the set 
is nearly uncorrected with its own shifted version while good 
cross-correlation means that any member of the sequences set 
is nearly uncorrelated with any other members at any shift. 

A commonly used metric of the goodness of the correlation 
is the Integrated Sidelobe Level (ISL). The ISL of a set of M 
sequences each of N (possibly complex) symbols that we will 
indicate with xj with j = 0, . . . ,N— 1 and p = 0, . . . ,M— 1 
is defined as 
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is is the auto-correlation of the sequence x^ p \ and 
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is the cross-correlation between the sequences x^ and x^. 

Good set of sequences are those having a low ISL value. 
Due to the strong interest in the design of sequences with 
low ISL value, many algorithms have been suggested for 
its minimization. Our purpose is to develop an analytical 
expression that may drive optimization in some particular 
difficult cases, most notably when the antipodal constrain 
(x'j = ±1) is imposed. 
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so that 
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A general method for the calculation of X jj,) (p) of any 
sequences of odd length is presented in [lj, [2J. This method 
hinges on generating functions and writes correlations as 
proper sums of their values on the unit circle in the complex 
plane. The method works well when we have analytical 
insights on the generating functions. 

Extending the ideas of JT|, in section [II] we devise a general 
method for the calculation of X jj,) ( ? ) in (|2]i of any pair of 
real sequences of odd length and thus, together with the 
result in JT], 0, the ISL for a set of sequences. In section 
III we use this method to obtain an asymptotic expression 
for the ISL value of a set formed by different rotations of 
Legendre sequences. Some minor results about an optimization 
procedure based on the latter expression are reported in (3), 
where we find the optimal rotations that minimize the ISL for 
any sequences length N. 

Throughout the paper we use the following asymptotic 
notation. 

We say that 

• two sequences and are asymptotically equivalent, 
a N ~ bpj iff 
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II. Calculation of the cross-correlation terms in 

THE ISL 

Let ao,ai, . . . ,fljv-i an( J £>0>^1>--->^V-1 be two real se- 
quences of length N, we want to obtain an expression for 

If we define the generating functions of the two sequences 

as 
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we have that 
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Now, set £: = e n J and note that for k,l = —N+l,. ..,N—l, 
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Hence, if we define 
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we can express X a & (i.e. the sum of squares of cross- 
correlations as in |2])) as 
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To compute 5" we use the Lagrange interpolation poly- 
nomials to calculate the values of Q a {— £/) from 2a( £ /<) f° r 
jf,t = 0, ...,N— 1. In this special case the data points (e^) 
coincide with the complex roots of unity and, for AT odd, the 
Lagrange base polynomials simply reduce to ^ ^ (4] p. 
89]. Then 
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in which we may exploit the fact that e* = 1 / Ej to write 
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Let us define now the innermost sum of |5) as 
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Depending on p,q,r,s, the rational function f Pt q <r ,s{z) can 
be transformed into a specific sum of simple rational parts. 
Each of these rational parts can be summed separately. This 
path is fully developed in JT] and we here exploit the results 
therein. 

In particular we have that 
A) for < p < N 
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By substituting Q into S" and developing the product 

\Q a {- £ j)Q* b {- £ j)\ 2 we get 
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Taking into account all the above cases we may write S" = 
M (oc + /3 + 7+5), where the terms a, /3, y, and 5 correspond 
to the contributions of the cases A, B, C and D respectively. 

For the cases included in A) we have that 
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for the cases in B) we have 
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for C) we have 
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Summarizing, we can write the sum of squares correspond- 
ing to cross-correlations terms of the ISL as 
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where the quantities a, j3, y, 5 are defined in |6]), 0, ((8), |9|. 

With the method presented above in conjunction with the 
method presented in [ 1 1, we can have an analytical expression 
for the ISL for any set of real sequences of odd length. The 
computation of the above equations seems to be hard at a first 
look, but in a number of cases, in particular for sequences 
from difference sets (2) may lead to significant results. 

In the following, we use this method to evaluate the 
asymptotic trend of the ISL of a set of sequences made up by 
different Rotations of a Legendre Sequence (RLS set) when 
N grows to infinity. 

III. Legendre Sequences 

The Legendre Sequence (LS) £0, . . . ,£n-i exists for any 
prime N and is defined as 



to 



1 if j is a square (mod N) 
— 1 if j is a nonsquare (mod N) 



A LS may be cyclically rotated t a positions to the left to 
obtain a Rotated Legendre Sequence (RLS) aj defined as 

a j — £j+t a (mod AO — ^j+f a N (mod N) 
With f a = ta/Ne [0,1]. 

The asymptotic value of X fla for the family of RLS was 
calculated in J5] and [j] noting that the asymptotic value of 
the modulus of the generating function of the LS (|G/'(£;)|) lS 
independent of j, yielding 
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We follow the same path as in U) but for the calculation of 
the cross-correlations terms of the ISL *K a b- 

To proceed, remember that the generating function of the 
LS is 



and for D) 
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{1 + ijVN if ./' 7^ and N =1 (mod 4) 
1 + UjVN if ./' ^ and N = 3 (mod 3) (11) 
1 if ./ = 

Moreover, if we denote by Q a (£j) the generating function 
of the RLS aj = l j+ta (mod N) , then 

Qa{£j) = £] ta Qt{£j) 

Assume now that the two sequences cij and bj are obtained 
by rotating £j by, respectively, t a and % positions to the left. 
We may compute S' as 



N-l 



N-l 



7=0 



S '=L £i ta QdejWj b QKej) =Z\Qdej 



N-l 



= 7,N 4 £ (X u (k)+X u (N-k)) 

° k=-N+l 
k^O 

c th + l _L c' h+2 _L c l ~' h _L c~' h -X- c'«+ 2 -X- -X- c 1- '" _l_ 

Note that X(j (k) + Xu (N — k) is the periodic correlation 
121 of the LS. Then, from and we know that 
\X(e(k) +Xa(N — k)\ < 3 for Legendre sequences. Then, using 
the fact that J^ 1 = 0{N 2 ) (see @ and (B) below 

and set t = 0), and using the triangle inequality we get that 
p = 0{N 6 ). 

• For the calculation of y in ([8j, following the same steps 
we did for j3 we have 
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from ( |TT| we know immediately that |<2^(£/)| 4 ^ N 2 , then 
S' ~N 3 , Let us now compute the asymptotic values of a, j3, 
7 and 5 in @, Q, ((8), <[§ for any pair of RLS. 
• For a in |6]) we have 
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and again we have that y= 0(N 6 ) 
< For 5 in |9} we have 
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Larger values of the summand are those for k close to 1, 
which make the denominator close to zero and numerator ~ cN 
for some constant c (for k close to N —I, the denominator 
becomes also close to zero but the numerator is 0(1)). 

Exploiting this and using the small angle approximation for 
the complex exponential, we may write 
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To continue, we recall the definition of the Dilogarithm 
function and its series expansion valid for \z\ < 1 
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Taking the real part of ( fT3] > and evaluating on the unit circle 
gives Q eq. (8.7)] 
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Exploiting ( fT4} and concentrating on the first period < 
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where = • (mod 1). 

Hence, since we know that 8 is real 
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where we have defined f a — ^ and f b = ^- Then, exploiting 
the symmetries of a quadratic form of a modulus function we 
have for < f a ,f b < 1 
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Fig. 1 . Plots of ISL for M = 2 as a function of fi and /*2 : (a) 3D- view, (b) 
iso-ISL lines 



Based on the above we are now interested in computing the 
asymptotic value of 
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Going back to our original problem for calculation of the 
ISL value of a set of M sequences Xj with j = 0,...,N—l 
and p = 0, . . . ,M— 1, where each jcW is made by a different 
rotation f p of a LS (RLS set), replacing ( fTO) and ( [TB) into ([3]) 
we finally have that 
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As an example, Figure [T] reports the 3D and contour plot of 
the right-hand side of (fT7|T for M = 2. Direct visual inspection 
of that Figure confirms that minima exists and can be easily 
identified. We offered a preliminary exploit of the result in 
where an optimization procedure was developed to find 
the optimal rotations that minimize the ISL for any sequences 
length N. 

As another example, in Figure [2] we plot the ISL for M = 4 
as a function of the sequence length N. In case a) the values of 
rotation are those that minimize the asymptotic ISL, while in 
case b) we use an arbitrary rotation. In both cases we can see 
that the trend of the plots is in agreement with the asymptotic 
value calculated. In the same figure in case c), we plot together 
both curves in a) and b) to show that the one that achieves the 
minimum asymptotic value of ISL, also achieves the minimum 
ISL value for sequences length greater than approximately 
20. For different choices of rotations and different number 
of sequences (M), the behavior is the same than presented. 

IV. Conclusion 

We apply a method based on generating functions, which 
has already been proposed for the calculation of the ISL 
of a sequence, to the calculation of the cross-correlation 
components of the ISL of a set of sequences. 

The apparent complexity of the resulting expressions can 
be tackled in the asymptotic conditions for sequences whose 
generating function has a relatively simple trend. 

Since this is the case of Legendre sequences, we are able 
to derive an analytical expression for the asymptotic ISL of 
sets of rotated Legendre sequences. 

Such an expression can be exploited to drive the optimiza- 
tion procedure needed to construct small-ISL sets of antipodal 
sequences with potential applications to communication and 
active sensing systems. 
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